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Abstract. Homogeneous isotropic models with two torsion functions built in the framework 
of the Poincare gauge theory of gravity based on general expression of gravitational La- 
grangian by certain restrictions on indefinite parameters are analyzed numerically. Special 
points of cosmological solutions at asymptotics and conditions of their stability in depen- 
dence of indefinite parameters are found. Procedure of numerical integration of the system 
of gravitational equations at asymptotics is considered. Numerical solution for accelerating 
Universe without dark energy and dark matter is obtained. It is shown that by certain 
restrictions on indefinite parameters obtained cosmological solutions are in agreement with 
SNe la observational data and Big Bang Nucleosynthesis predictions. Statefinder diagnostics 
is discussed in order to compare considered cosmological model with other models. 
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1 Introduction 

One of the most principal recent achievements of observational cosmology is the discovery 
of the acceleration of cosmological expansion at present epoch [1, 2]. In order to explain 
this observable accelerating cosmological expansion in the framework of General Relativity 
Theory (GR), the notion of dark energy (or quintessence) was introduced in cosmology. 
According to obtained estimations, approximately 70% of energy in our Universe is related 
with some hypothetical kind of gravitating matter with negative pressure — "dark energy" 
— of unknown nature. Previously a number of investigations devoted to dark energy problem 
(DEP) were carried out (see reviews [3, 4]). 

According to widely known opinion, the dark energy is associated with cosmological 
term, which is related in the framework of standard ACDM-model to the vacuum energy 
density of quantized matter fields. In connection with this the following question appears: 
why the value of cosmological term is very small and close to average energy density in 
the Universe at present epoch (see for example [5]). Other treatment to solve the DEP is 
connected with modification of gravitation theory. 

One of such solutions can be obtained in the framework of gravitation theory in the 
Riemann-Cartan spacetime C/4 - Poincare gauge theory of gravity (PGTG) [6-10] . It should 
be noted that the PGTG is a natural and in certain sense necessary generalization of metric 
theory of gravitation by including the Lorentz group to the gauge group which corresponds 
to gravitational interaction. The PGTG leads to the change of gravitational interaction 
in comparison with GR and Newton's theory of gravity at cosmological scale, which are 
provoked by more complicated structure of physical spacetime, namely by spacetime torsion 
[11, 12]. 



Explicit form of gravitational equations of PGTG and their physical consequences de- 
pend essentially on the structure of gravitational Lagrangian vCg, which is built by means of 
invariants of gravitational gauge field strengths — the curvature and torsion tensors. The 
most simple PGTG is Einstein-Cartan theory based on the gravitational Lagrangian in the 
form of scalar curvature of spacetime C/4 [13]. In the frame of Einstein-Cartan theory the 
DEP was discussed in [14], where some phenomenological description of spinning matter was 
used. In connection with this it should be noted that in the frame of Einstein-Cartan theory 
the torsion is connected with spin momentum by linear algebraic relation and vanishes in the 
case of spinless matter. Such situation seems unnatural by taking into account that the tor- 
sion tensor is gravitational gauge strength corresponding to transformations of translations, 
which are connected directly with energy-momentum tensor in the frame of Noether formal- 
ism. The situation comes to normal by including to £g terms quadratic in the curvature and 
torsion tensors. 

The most general form of the gravitational Lagrangian Cg (without using Levi-Civita 
symbol) includes the linear in the scalar curvature term as well as 9 quadratic terms (6 invari- 
ants of the curvature tensor and 3 invariants of the torsion tensor with indefinite parameters) . 
The structure of the gravitational equations and physical consequences of isotropic cosmology 
in the frame of PGTG, in particular, the situation concerning the DEP depend essentially 
on restrictions on indefinite parameters. The PGTG can be divided into diff'erent sectors in 
dependence on the number of nonvanishing componets of the torsion tensor and the order of 
the differential equations determining their behaviour and behaviour of the scale factor that 
is connected with restrictions on indefinite parameters of £g. 

In general case the torsion tensor for homogeneous isotropic models (HIM) is described 
by two functions of time Si and S2 determining trace and pseudotrace of the torsion tensor 
respectively. At the first time the most simple HIM with the only nonvanishing torsion 
function Si were built and investigated in [15]; it was shown that by certain restrictions 
on equation of state of gravitating matter at extreme conditions (extremely high energy 
densities and pressures) in the beginning of cosmological expansion all cosmological solutions 
are regular with respect to metrics, Hubble parameter, its time derivative and energy density 
by virtue of gravitational repulsion effect at extreme conditions (see [16]). The regular Big 
Bang scenario with infiationary stage in the beginning of cosmological expansion based on 
such HIM was built and analyzed in [17, 18]. 

Other sector of PGTG is so-called dynamical scalar torsion sector considered in [19, 20]. 
HIM biult in this sector demonstrate the oscillating behaviour of the Hubble parameter, and 
it is possible to obtain good correspondence with SNe la observational data. 

By investigation the sector of PGTG with two torsion functions it was shown that the 
PGTG allows to explain the acceleration of cosmological expansion at present epoch without 
using the notion of dark energy [21].^ This result is due to the fact, that the cosmological 
equations for HIM at asymptotics take the form of Friedmann cosmological equations of 
GR with effective cosmological constant induced by spacetime torsion if certain restrictions 
on indefinite parameters of gravitational Lagrangian are imposed. As it was shown in [25], 
isotropic cosmology based on HIM with two torsion functions offers opportunities to solve 
also the problem of dark matter. The analysis of regular inflationary cosmological models 



^At the first time equations for HIM with two torsion functions were deduced in [22]. These equations 
were considered in [23] with the purpose to obtain their solutions; however, so called "modified double duality 
ansatz" used in [23] by obtaining solutions with non-vanishing torsion function 5*2 is not applicable in this 
case even for the vacuum (see [24]) and its application leads to incorrect solutions. 



built on the base of HIM with two torsion functions by certain restrictions on indefinite 
parameters in gravitational equations for such models was fulfilled in [26]. The present paper 
is devoted to numerical analysis of HIM with two torsion functions of accelerating Universe 
at asymptotics when energy densities are sufficiently small. 

2 Cosmological equations for homogeneous isotropic models 

In this Section we briefly repeat the derivation of the cosmological equations for HIM with 
two torsion functions (see [21, 26]). 

In the framework of PGTG the role of gravitational field variables play the tetrad /i'^ 
and the Lorentz connection A^^fj^; corresponding field strengths are the torsion tensor 5*^,^ 
and the curvature tensor F ^j^y defined as 

S\v = d[y h'^] - hk[fj,A' „] , 

F'\, = 2d^^A\^+2A'\^A\i\,^, 

where holonomic and anholonomic space-time coordinates are denoted by means of greek and 
latin indices respectively. 

We will consider the PGTG based on gravitational Lagrangian given in the following 
general form 



£g = /o F + F"^^'^ (/i F^p^, + /2 F^^p, + /3 F^,^p) + F^^ [U F^u + h F, 

+h F^ + 5"^^ (ai S^^, + as S,^o) + ^3 5"^,5^' 



.5/^ (2.1) 

where F^y = F^^au, F = F^^, fi (i = 1, 2, . . . , 6), a^ {k = 1, 2, 3) are indefinite parameters, 
/o = (167rG)~^, G is Newton's gravitational constant (the velocity of light in the vacuum 
is equal to 1). Gravitational equations of PGTG obtained from the action integral / = 
J {Cg + Cai) hd'^x, where h = det (/i*^) and C„i is the Lagrangian of gravitating matter, 
contain the system of 16+24 equations corresponding to gravitational variables /i*^ and A fj_. 
The sources of gravitational field in PGTG are the energy-momentum and spin tensors. In 
present paper we will consider perfect fluid with energy density p, pressure p = p{p) and 
vanishing spin tensor as a source of gravitational field. 

High spatial symmetry of HIM allows to describe these models by three functions of 
time t: the scale factor of Robertson- Walker metrics R and two torsion functions Si and 5*2 
determining the curvature functions ^4^ {k = 1,2,3,4) as following 

Ai = H + H'^ - 2HSi - 25*1, 

A2 = ^ + {H-2Sif-Sl 

A3 = 2{H-2Si)S2, 
Ai = 5*2 + HS2, 

where H = R/R is the Hubble parameter and a dot denotes the differentiation with respect 
to time. The system of gravitational equations of PGTG for HIM in considered case takes 



the following form: 



, (i7 _ s,) Si - 2bSl - 2/0^2 + 4/ {Al - Al) + 2(/2 (^1 - ^1) = -^, (2.2) 



where 



' -Al) 

2 yl2^ 



l2 a2\ 



Si + 2HSi -Sf)- 2bSl - 2/o (2Ai + A2) - 4/ (^^ _ ^2^ _ 2^3 (^f - Al) = p, (2.3) 



/ 



92 



a 



Al + 2/7 {Al - A2) + 4SiA2\ + 92^2^3 - qiS2A^ + (^/o + g j 5i = 0, (2.4) 
^4 + 2H (^4 - A3) + ASiA-i] - 4/ S2A2 - 2qiS2Ai - (/o - 6) ^2 = 0, (2.5) 



2ai + a2 + 3a3, 

/2 



b = a2 — ai, 



/ = /i + y + /3 + /4 + /5 + 3/6 , 
9l=/2-2/3 + /4 + /5 + 6/6, 92 = 2/1-/2. 

The system of gravitational equations (2.2)-(2.5) allows to obtain the cosmological equations 

generalizing Friedmann cosmological equations of GR and equations for the torsion functions 

5i and ^2. 

To exclude higher derivatives of the scale factor R from cosmological equations the 

following restriction on indefinite parameters a^ was imposed: a = (see [15, 24]). Obtained 

gravitational equations can be simplified if an additional restriction on parameters fk is 

imposed 2f = qi + 52 [21]. Then cosmological equations and equations for torsion functions 

contain three indefinite parameters: parameter a = o -4- > with inverse dimension of energy 

■^ /o 
density, parameter b with dimension of parameter /o and dimensionless parameter e = y-. 

The particle content of the PGTG with these restrictions on indefinite parameters of the 

gravitational Lagrangian (2.1) was discussed in ref. [26]. 

For further analysis, we transform cosmological equations to dimensionless form by 

introducing dimensionless units for all variables and parameter b entering these equations 

and denoted by means of tilde: 



t^i = t/y/6f^, R 

P^P=^ap, p- 

S'1,2 -^ <Si,2 = 'S'i,2\/6/oa, b - 



> R = R/y/6ha, 
p = ap, 
b = b/fo, 



(2.6) 



H -^ H 



where dimensionless Hubble parameter H is defined by usual way H = Rr^^. As result 
cosmological equations ((22)-(23) in Ref. [21]) take the following dimensionless form, where 
the differentiation with respect to dimensionless time t is denoted by means of the prime and 
the sign of "" is omitted below in this Section and Sections 3 and 4: 



k_ 
R? 



+ iH-2Si] 



1 
Z 



p + {Z-b)Si + ^{p-3p-2bSl)^ 



£ 

2Z 



{HS2 + s'y + ^ 



k 

i?2 



S^ ] S? 



(2.7) 



H' + H"^ - 2HSi 



-^{p-^P 



-'''^ = -2Z 

2bSl) Si + — 
2; 2 2Z 



P + ^P- ^\P 



2,p-2bSlY 
{HS2 + S'2)' + 4(^-^-Sl^S: 



(2. 



(Z = 1 + p - 3p - 2 (6 + e) 5|) . 

In considering case of HIM filled with usual gravitating matter with equation of state in the 
form p = p{p) the torsion function Si in dimensionless form appearing in (2.7)-(2.8) is 



S. = -J^<H 



(p + p)(3|-l)+2.5? 



^+2eSl - I {2b -e) 828'^'^ (2.9) 



and dimensionless torsion function 52 satisfies the following differential equation of the second 
order: 



e [8'^ + 3H8'2 + 3H'82 - 4 (5^ - 3H8i + ASf) 82] 
-2 {p-3p- 265|) ^2 - 2 (1 - 6) S2 = . 



(2.10) 



The conservation law for gravitating matter in dimensionless units has the usual form 

p' + 3H{p + p) = 0. (2.11) 

3 Critical points analysis 

The system of equations (2.8) - (2.10) together with conservation law (2.11) completely 
determine the dynamics of HIM, if the equation of state of matter is given. For further 
analysis we will consider flat model {k = 0) filled with matter with barotropic equation of 
state p = wp {w = const). The aforementioned system of equations can be represented in 
the form of four first order differential equations for H, 82, U = 82 and p: 



MoY' = F, 



where the matrix Mq is 



Ma 



and 



/ 1 


2^s^ 



-2i| -2ff - 
1 


-2f^ 



3e82 - 4e^52 

V 


-4ei|52 -4.f^52 -4.^52 
1 


Y = 


(H\ 

82 

U 

\P ) 


J 


F = 


/:FiiH,82,U,p)\ 
:F2iH,82,U,p) 
T3{H,82,U,p) 

\MH,82,U,p)J 


; 



(3.1) 



(3.2) 



J^i 



T2 

^4 



i[(l-3u;)/,-26S|]'} 
-U(x-3w)p-2{h-\)8''2\8l^^{H82^U)\ 



-H^^2H8i- — \{\^3w)p 



(3.3) 

[/, (3.4) 

-e [3HU + 4 {3H8i - A8l) 82] + 2 [(1 - 3w) p - 2b8l] ^2 + 2 (1 - 6) ^2, (3.5) 
-3{l + w)pH. (3.6) 



The function Si can be written as 
Si 



_3_ 



Ih [(1 + w) {3w -l)p + 2eS|] - ^ (2fe - e) S2u\ . (3.7) 



Critical points Pi = Pi{Hc, S2c, Uc, Pc) of the first order system of differential equations 
(3.1) can be obtained by setting H' , S2, S[, p' to zero [27, 28], i.e. by solving the following 
system of equations: 

J^,{H,S2,U,p)=0, (i = l,...,4). (3.8) 

In the case of flat model (k = 0), solutions of (3.8) have to satisfy (2.7) with A; = 0. Eq. (3.4) 
leads to [4 = 0. 

Obviously, the point Pq with vanishing values of H^ S2C, Pc satisfies (3.8). Analogously 
to GR this point is the point of complicated equilibrium. To analyze the stability of other 
critical points P{Hc, S2c,0, Pc) satisfying (3.8) it is necessary to build linearized form of the 
system (3.1). Near the critical point the variables can be written in the form H = Hc + Aif , 
S2 = S2C + A52, U = AU, p = Pc + A/3 and the linearization of the system (3.1) takes the 
following relation 

AY' = Mq^Mi ay, (3.9) 

where the components of the matrix Mi are given by 



Ml 



i,*i 



dYj 



Stability of the point P is determined by the eigenvalues Aj of the matrix M^ Mi [27, 28] . 
Characteristic equation det (Mi — AMq) = leads to quartic expression with respect to A, 
which can be written as 

0. (3.10) 



A'^ + CiA^ + C2A^ + C3A + C4 



If the real parts of all Aj is negative, then the critical point P is stable and the gravitational 
equations (2.8) - (2.11) can have asymptotics to this point H — )• He, S2 — t- S2c, S2 — )■ 0, 
/9 —7- /9c at t —7- +00. 

According to the Routh-Hurwitz theorem all Aj will have negative real parts if the main 

minors of the matrix 

/ ci 1 \ 

C3 C2 Ci 1 
C4 C3 C2 

V C4 / 



(3.11) 



are positive [27], i.e. 

ci > 0, C1C2 - C3 > 0, 



C1C2C3 - C1C4 - c| > 



and 



C4 > 0. 



(3.12) 



The equation J^4(Hc,S2c,0, pc) = gives two kinds of critical points: with vanishing 
Hubble parameter He = or with vanishing energy density pc = 0. Let us consider in details 
both of them. 



3.1 Critical points with vanishing Hubble parameter 

If He = the system (3.8) is reduced to the system of two algebraic equations 

^ I (1 + 3u;) p - ^ [(1 - 3w) p - 2bSl]^\ + e [{1 - 3w) p - 2(1 - 1) 5|] 5| = 0,(3.13) 

[l-b+{l-3w)p-2bS^]S2 = 0. (3.14) 

Generally speaking, except trivial solution S2c = 0, pc = this system admits non-zero 
solutions for S2c and pc- According to (3.14) there are the following possibilities: 52 = and 
5*2 7^ 0. In the first case we obtain the trivial solution with pc = 0. If S2c 7^ 0, from (3.14) it 
follows that (1 — 3w)p = 6 — 1 + 2651. As result we have 

and Z = b-2eSl (3.15) 











Si = 





Substitution 


(3.15) 


into 


(2.7) 


gives 










k 


/ 

1 + 
V 


2e 
~Z 



1 



p+\{b-\f 



0, (3.16) 



and this equation does not have non-trivial solution in physical space (p > 0) for flat models 
(fc = 0). 

3.2 Critical points with vanishing energy density 

If /9c = the system (3.8) is reduced to the system of two algebraic equations 

-U^ + 2iISr + |5| + ^ii^5| + ^R'Sl = 0, (3.17) 

[2e (3F5i - 45?) + 265| - (1 - b)\ S2 = 0, (3.18) 

where the functions 5i and Z can be represented in the following form 

Si = -^HS^ and Z = 1 - 2 (6 + e) 5|. (3.19) 

ZZj 

Neglecting the case 52 = 0, it is possible to obtain from (3.17) - (3.18) equation for H 
in closed form. To do this, the system of equations (3.17)-(3.18) can be rewritten in the 
following form 

2(1- 265|) H^ + eH^Si - 2 [6^ + 2e {b - 1)] 5| = 0, (3.20) 

9e[l-2{b + e- e^)] H^S^ - 2bS^Z^ + {1 - b) Z^ = 0. (3.21) 

Polynomials in the left-hand-side of equations (3.20)-(3.21) generates ideal in the polynomial 
ring in the variables H and 52 [29]. There are different ways to choose the basis in the 
ideal of polynomials. One of them is a Grobner basis with lexicographic ordering H ~< 82- 
Using computer algebra system Wolfram Mathematica the first element of Grobner basis of 
aforementioned polynomial ring takes the form 

(-9e(5 + 2e)H^ + b^ + 2be - 2e) {2be^ [Sb^ - be{9e + 2) + 86^] H^ 

-e [46^(e - 2) + b^e{l7e + 22) - b'^e {32e^ + 41e + 34) + Sbe^i^e + 3) - 32e^] H^ 

-2 [b^{£ + 4) + 6^(5 - 6e)e - 2b^e [Se^ - 4e + 3) + 2b^e^{2A£ + 1) 

-46e2(i2e + 1) + 16e^] H'^ + 2(6 - 1)^ (6^ + 26e - 2ef\ = 0. (3.22) 



Repeating this procedure with lexicographic ordering S2 ~< H we have the first element of 
Grobner including only 5*2 

[2(6 + 2e)Si -I] {2b [8b^ - be{9e + 2) + Se^] sl + [8b^ - 2b^{e + 12) 

+46e(2e + 1) - Se^] 5| - [86^ + b{e + 12) + e] 5| + 2(6 - 1)} =0. (3.23) 

According to general mathematical theorems [29] roots of the system of equations (3.20)- 
(3.21) turn (3.22) and (3.23) into true. 

Analytic analysis of stable points determined by the system (3.17)-(3.18) is possible 
approximately only if 1 — 6 — )• +0 and e — )■ 0. In other cases it is necessary to use numerical 
methods. 

3.2.1 Approximate analysis in the case < 1 — 6 ^ 1 

System of equations (3.17) - (3.18) admits simple approximate solution if < 1 — 6 ^ 1. 
This solution was initially obtained in ref. [21] and after transformation to dimensionless 
form (2.6) it reads: 

This approximation is valid up to cubic term in 1 — 6 and linear term in e. 

In the case < 1 — 6 <C 1 the stability of the critical point P2 ~ ( iFtt-, y ^is") Oi ) can 
be analyzed analytically. In this case the inequalities (3.12) leads to: 

e > 0, w> -I. (3.25) 

3.2.2 Numerical analysis of stability 

As an exact analytic expression for solution of the system (3.17)-(3.18) does not exist in 
general case, it is necessary to use numerical methods to analyze stability of the critical 
points. The procedure of the numerical analysis of the stability points is following. 

1. For given value of e, the system (3.17)-(3.18) is solved numerically for the set of values 
6. 

2. For every real solution of the system (3.17)-(3.18) at given values of e and 6 character- 
istic equation det (Mi — AMq) = has to be solved with respect to A. 

3. The real parts of obtained Aj have to be tested for negativity. 

For example, the results of this procedure for e = 0.03 are given in figure 1. The calculation 
are performed for 6 varying from 0.01 to 1.2 with a step A6 = 0.05. In the left panel of 
figure 1 the curves determined by (3.22) are imposed. In the right panel an analogous curves 
for 52 determined by (3.23) are imposed. 

From figure 1 it is possible to see, that there is minimal value of 6 assuming nontrivial 
solution of eqs. (3.17)-(3.18). This value can be found by setting H to zero in (3.22). As 
result we have the following restriction on 6 



6 > -e + ^/e{2 + e). (3.26) 





Figure 1. He and 820 as function of 5 (e = 0.03). Stable critical points are marked by (green) circle. 
Unstable critical points and points of complicated equilibrium are not shown in this figure. Solid lines 
in the left and right panel of figure correspond to solutions of eqs. (3.22) and (3.23) respectively. 



4 Numerical integration of the system of gravitational equations 

In this paper we will analyze the late time behaviour of the solution of the system (2.8)- 
(2.11). To make comparison with GR we will perform numerical integration of the system 
of the gravitational equations for dust matter (tt; = 0) . To simulate late-time behaviour the 
initial conditions will be taken at to = according to the following procedure. 



1. For given value of e and 6 (e > 0, — e + \Je (2 + e) < 6 < 1) algebraic system (3.17)- 
(3.18) is solved numerically and all critical points Pj(iJc, 5'2c,0, 0) are found. Only real 
solutions are considered. 

2. For every critical point, the stability analysis is carried out according to the previous 
subsection and stable point with minimal positive He and positive S2C is selected.^ 



3. The torsion function 5*2 and the Hubble parameter H can be represented in the form 

(4.1) 
(4.2) 



sl = sic + y2 p, 



with some coefficients yi and ^2-'^ As the stable point is selected, then p tends to zero 
at t — 7- +00. Keeping linear terms in p the conservation law (2.11) can written as 



-SHcp. 



(4.3) 



Substitution of (4.1)-(4.3) into (2.8)-(2.10) together with keeping terms linear in p 
gives two algebraic equations for determination of yi and y2- Numerical solution of 
these algebraic equations for given e, b, He and 52c gives yi and 7/2- 

^These values of He and ^c correspond to the vacuum as de Sitter spacetime with torsion [24]. 
^Further in this paper we will refer representation (4.1)-(4.2) of H and S2 as late-time approximation. 




Figure 2. Late-time behaviour of Hubble parameter and S2 torsion function. 





Figure 3. Late-time behaviour of the deceleration parameter and energy-density. 



4. Positivity of obtained values of yi and 2/2 is considered as applicability of the late time 
approximation (4.1)~(4.2) and successful choice of stable critical point made in step 2 
of current procedure. Further steps are performed only if yi > and ^2 > 0. 



5. Initial condition for po = p(io) is taken from the following equation 

H^to) _ H^ + yipo 1 



m 



m 



^A 



(4.4) 



as result we have Hq = H{to) = \/H^ + yi po and ^20 = »S'2(to) = v^lc + ysPo- Here 
Qa is an additional free parameter that specifies initial conditions. 

6. Initial condition for 82^ = »S'2(to) is obtained from (2.7) taking into account A; = 0. The 
minimal in modulus value of iS'20 is taken as initial value. 

7. For this choice of the parameters e, h and initial conditions /jq, -f^O) 'S'20 and 6*20 the 
system of differential equations (2.8)-(2.11) is integrated numerically. 

As an example let us consider the numerical solution at the following parameters and 
initial conditions e = 0.00019, h = 0.14, Ho = 1.3417, S20 = 1.82574, S'^q = -0.288424, 



10 



po = 0.0730922. This choice of the initial conditions gives H'^ (to) / H^ (oo) = I/^a = 1/0.72. 
Figures 2-3 show the behaviour of Hubble parameter H, torsion function 5*2, deceleration 
parameter 

q = -R"R/R'^ (4.5) 

and energy-density of dust matter p. Another choice of the initial conditions leads generally 
speaking to cosmological solutions with another behaviour of the Hubble parameter and 
torsion function ^2 because of their oscillating character. 

5 Comparison with observational data 

In this Section we will analyze what restrictions on indefinite parameters leads to solutions 
corresponding to observational data. It should be noted that the numerical solution of 
gravitational equations allows to obtain time dependence of Hubble parameter H{t) and scale 
factor R{t) for given values of parameters of gravitational Lagrangian and initial conditions. 
As result dimensionless luminosity distance di, can be obtained as a function of redshift z 

, = :^_i = 5M_i (5.1) 

R{t) R{t) ^ ' 

in the following form [30, 31]: 

Jl = R{to){l + z) r jL = ^i + ,) r Jl_ (5.2) 

Ji R{t) Jo H{z) 

The predicted distance modulus fi = m — M {m and M are apparent and absolute 
magnitude respectively) can be written as function of dimensional luminosity distance c^l = 
di,y/6foa in the following form 

/i = m-M = 25 + 5 logio^L, (5.3) 

where c/l is given in megaparsecs. 

5.1 Matching the late time approximation 

At first we will compare late time approximation (4.1) of H{t) with supernovae type la 
(SNe lA) observation data and predictions of the standard Big Bang Nucleosynthesis theory 
(BBN). Dependence of the energy density of the dust matter p{t) ~ R^it) allows to write 
Hubble parameter in the late time approximation (4.1) as function of the redshift z. As 
result the corresponding expression for the predicted distance modulus /i takes the form 



p{z) = 25 + 5 log 
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Jo 



dz 



where 

/io = 5 logio 



A/6/oa 



+ /io, (5.4) 



Hq ■ IMpc 



11 



5.1.1 Comparison with SNe la observational data 

We will start from comparison with supernovae observations using the Union2 compilation 
of 557 SNe la data [32] (see also [33]) and minimize 

XSN — 2^ 2 ' v^-''/ 

where dj is the distance modulus errors. 

Best fit parameters Q,a = 0.73 and /.to = 18.16 gives Xsn ~ 542.683 and Xsn/'^^^ ~ 0.974 
(dof — degree of freedom). The value of /io = 18.16 corresponds to 70km/(sec • Mpc) for 
Hubble constant at present epoch as in ACDM-model of GR. Indeed, late time approximation 
(4.1) of H{t) is similar to Friedmann equation of GR, but differs only by effective gravitational 
constant Gcs determined as Gcs = ViG. 

If the matter content includes baryonic and dark matter with relative contributions fie 
and r^DM to the total energy density, than it is easy to show [25] that 

yi = 7^ — —pp- (5-6) 

In particular, as matter candidate for cold dark matter is not found yet, it is possible to fit 
SNe la observational data in the discussed model without using dark matter (0dm = 0). 

5.1.2 Comparison with SNe la + BBN data 

Calculations in the frame of standard BBN theory predicts ujb = ^b^^ = 0.0212 it 0.0010 
for baryon mass density [34], where h is the Hubble constant at present epoch in units of 
100km/(secMpc). For h = 0.7 obtained earlier the corresponding Ob = 0.043. Assuming 
that the rate of light element production during Big Bang Nucleosynthesis does not depend 
on the presence of the torsion^ and the dynamics of H{t) from BBN epoch to present epoch 
is well approximated by (4.1), equation (5.6) gives yi = 6.3 for Q\ = 0.73 and 0dm = 0. As 
yi is completely determined by parameters b and e, it seems impossible to determine b and 
e simultaneously using late time approximation (4.1). 

5.2 Matching general case 

Besides approximation (4.1) there is another way to obtain the best fit parameters of consid- 
ered theory. Namely, for given values of e, b and specified initial conditions, the procedure of 
numerical integration of exact system of differential equations (2.8)-(2.11) allows to obtain 
solution for H(t), p(t) and predicted distance modulus as a function of redshift // = iJ.{z). Ob- 
tained distance modulus - redshift dependence n{z) allows to calculate joined x^ for Union2 
data set and BBN predictions 

X^ = Xsn + xIbn> (5-7) 

where 

-.2 



2 _ (ObIO-'^o/^c/ (100 km/sec) - cjb) 

XbBN — ~2 

VB 



This assumption does not contradict to results obtained in [35]. 
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c is the velocity of light and for computational purposes the functions fj,(z) and Ob are written 
in the following form 



IJ.{z) = 25 + 5 log 



10 



(1 + ^)/ 
Jo 



dz 



H{z) 



+ fJ-o, (5. 



_ 1-0a 

i'B — "DM- 

yi 

According to [34] we will use ujb = 0.0212 and cjyg = 0.0010. 

In general case x^ depends on parameters a, b, e, initial conditions po, Hq, S20, 820 {^a) 
and cold dark matter contribution J^dm to the total energy density. The task of minimization 
of total x^ implies the minimization with respect to all parameters and initial conditions. To 
simplify this problem we will restrict the task by setting initial conditions in dependence on 
the values of parameters as was discussed in Section 4. 

Following the previous subsection at first we will use the Union2 compilation set only. 
Considering the grid in the parameter region 0.00009 < e < 0.0002, 0.12 < b < 0.2, 0.66 < 
Oa < 0.74 (Ae = 10~^, Ab = 0.01, A^a = 0.01) and minimizing Xsf^ we find the minimum 
of xIn at e w 0.00019 and b « 0.12, Oa ~ 0.72 and flo = 19.0074 (xIn = 542.8, XsN/^iof = 
0.975). Obtained values of parameters correspond to 70.0km/(sec-Mpc) for Hubble constant 
at present epoch, ujb = 0.0173 and a~^ = 4.22 • 10"^*^ g/cm'^. Obtained value of wb seems to 
be small and contradicts data on D- and ^ He- abundance, but close to data on ^He [36]. 

Calculation x^ defined by (5.7) in the same grid in parameter region 0.00009 < e < 
0.0002, 0.12 < b < 0.2, 0.66 < f^A ^ 0.74 we find approximate best fit parameters for this 
model: e ^ 0.00019, b ^ 0.14, ^a ^ 0.72 and /2o = 18.799 (x^ = 542.9, xVdof = 0.983). 
Obtained values of parameters correspond to 69.9km/(sec • Mpc) for Hubble constant at 
present epoch, a~^ = 5.11 • 10~^''g/cm^ and ojb = 0.0198 which is in accordance with data 
on ^ He- abundance, and lies in 2a interval for D-abundance. Solution presented in Figures 2- 
4 corresponds to this set of indefinite parameters. Comparison of the dependence of the 
distance modulus /x as a function of redshift z for obtained numerical solution with that in 
ACDM-model and Union2 observation data is presented in Figure 4. 

6 Statefinder diagnostics 

It was shown in a number of papers [37, 38] that so-called statefinder diagnostics proposed by 
Sahni, Saini, Starobinsky and Alam [39] allows to effectively discriminate between different 
models of dark matter and dark energy using the data from future SNAP-type satellite 
missions [40]. 

Statefinder diagnostics was applied to cosmology based on dynamic scalar torsion sector 
of PGTG [41] and it was found that some characteristics of the evolution of statefinder 
parameters 

r = R"'/RH^ and s = (f - l)/3(g - 1/2) (6.1) 

can be distinguished from that for other cosmological models.^ The evolutionary trajectories 
of the statefinder pair (f, s) for numerical solution obtained in Section 4 are shown in the 



The statefinder parameter r is also known as jerk j [42]. 
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Figure 4. Comparison with SNc la observational data. Solid line corresponds to obtained numerical 
solution and dashed line corresponds to ACDM-model. 





Figure 5. Evolutionary trajectories of statefinder pair 

Figure 5 in the {q, f) and (s, f) planes. This Figure demonstrates that statefinder diagnostics 
allows to distinguish cosmology in considered sector of PGTG from cosmology based on 
scalar torsion sector of PGTG [41] and other cosmological models, that allows in principle to 
differentiate the considered cosmological models from others using future SNAP-type satellite 
missions. It is necessary to note complicated and oscillatory behaviour of statefinder pair 
which is result of oscillating behaviour of the deceleration parameter (see Figure 3). As 
result the successful comparison observation data from planned SNAP-missions may require 
particular procedure for observation data processing. 

One of the possible procedure may consist in calculation of the statefinder pair (f, s) 
based on data from specially selected intervals (zj, Zj+i). For example, averaging of the 
statefinder pair for the obtained numerical solution in the range 0.069 < z < 0.26 gives 
(f) = 0.88 and (s) = 0.057 which is close to that for ACDM-model (f, s)acdm = (1, 0), 
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Figure 6. Statcfindcr pair as function of rcdshift z. 

but averaging over intervals 0.069 < z < 0.16 and 0.16 < z < 0.26 gives ((f) , (s)) equal to 
(2.05, —0.36) and (—0.15, 0.43) correspondingly. Thus, averaging over twice smaller interval 
demonstrates different values and oscillations near the point (1, 0)). This feature in the 
behaviour of the state finder pair can be used as a crucial test of the considered model in the 
planned SNAP-type satellite missions, but comprehensive analysis of possible tests of such 
type including the procedure of Zi determination is out of this paper. 

Conclusion 

As follows from our analysis, homogeneous isotropic models built in the framework of the 
Poincare gauge theory of gravity and filled by ideal fluid can have stable solutions with de 
Sitter asymptotics, if certain restrictions on indefinite parameters of the gravitational La- 
grangian are imposed. Obtained model demonstrates accelerated expansion at the late time 
approximation and does not include dark energy for which physical nature is still unknown. 
Contribution of dark matter to the matter energy density is a free parameter of considered 
model and it can be either vanishing or non-vanishing. In the case of vanishing dark matter 
only baryonic matter with dust equation of state contributes to the total energy density in 
considered model. 

Correspondence with SNe la observation data and BBN predictions are analyzed, where 
restricted set of indefinite parameters was considered and special procedure for determination 
of an initial conditions was used. Using this procedure best fit values of indefinite parameters 
and initial conditions are found. Obtained numerical solution was shown to be in a good 
correspondence with ACDM-model as well as SNe la observation data and in an accordance 
with data on "^He- and D-abundance. 

It was shown, that the trajectories of the statefinder pair for obtained solution demon- 
strate behaviour different from that in dynamic scalar torsion sector of PGTG and other 
cosmological models, that allows in principle to discriminate the considered cosmological 
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models from others. Special feature in the behaviour of this trajectories is noticed allowing 
to test considered model using data from future SNAP-type satellite missions. 
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